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In previous work MajciaEi0 and McLaughli MB computed explicit expressions for the 2Nth moments of a 
passive scalar advected by a linear shear flow in the form of an integral over . In this paper we first compute 
the asymptotics of these moments for large moment number. We are able to use this information about the large 
N behavior of the moments, along with some basic facts about entire functions of finite order, to compute the 
asymptotics of the tails of the probability distribution function. We find that the probability distribution has 
Gaussian tails when the energy is concentrated in the largest scales. As the initial energy is moved to smaller 
and smaller scales we find that the tails of the distribution grow longer, and the distribution moves smoothly 
from Gaussian through exponential and "stretched exponential" . We also show that the derivatives of the scalar 
are increasingly intermittent, in agreement with experimental observations, and relate the exponents of the scalar 

derivative to the exponents of the scalar. 
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I. BACKGROUND 

It is a well documented experimental fact that, while the statistics of the velocity field in a turbulent 
flow are roughly Gaussian, the statistics of other quantities like the pressure, derivatives of velocity and 
a passively advected scalar are generally far from Gaussian. [[7|J^ JT5| , |l^ , p3| , pO| For example Castaing, et. 
al. ||] observed in experiments in a Rayleigh-Benard convection cell that for Rayleigh number Ra < 10^ 
the distribution of temperature appeared to be roughly Gaussian, while for larger Rayleigh numbers, 
Ra > 10^, the temperature distribution appeared to be closer to exponential. In related work Ching | ]l5t 
studied the probability distribution functions (pdfs) for temperature differences at different scales, again 
in a Rayleigh-Benard cell, and found that the pdfs over a wide range of scales were well approximated 
by a 'stretched exponential' distributions of the form 
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P{T) = e-^l^l". 

At the smallest scales the observed value of the exponent was /3 ~ .5, while at the largest scales the 
observed exponent was roughly /3 « 1.7. Kailasnath, Sreenivasan and Stolovitky measured the 
pdfs of velocity differences in the atmosphere for a wide range of separation scales. They found similar 
distributions to the ones found by Ching, with exponents ranging from /3 w .5 for separation distances 
in the dissipative range to /3 w 2 on the integral scale. Finally Thoroddsen and Van Atta ||5^ studied 
thermally stratified turbulence in a wind tunnel and found the probability distributions of the density to 
be roughly Gaussian, while the distributions of the density gradients were exponential. 

A complete understanding of such intcrmittency lies at the heart of understanding fluid turbulence, and 
would certainly require a detailed understanding of the creation of small scale fluid structures involving 
both patchy regions of strong vorticity and intense gradients p^Ji^ . An alternative starting point is to 
assume the statistics of the flow arc known a priori and to determine how these statistics are manifest in a 
passively evolving quantity. This question of inherited statistics is significantly easier than the derivation 
of a complete theory for fluid turbulence, though still retains many inherent difficulties such as problems 
of closure. 

Motivated by the Chicago experiments of the late 80's and earlier work |^,^ 36 4^, there has been a 



tremendous effort towards understand the origin of the intermittent temperature probability distribution 
function in passive scalar models with prescribed (usually Gaussian) velocity statistics. For a very 
complete review of the subject of turbulent diffusion, including a full discussion of scalar intcrmittency, 
see the recent survey article of Majda and Kramer |3|]. Most of the work on the scalar statistics has 
either been directed at understanding the anomalous scaling of temperature structure functions, or at 
understanding the shape of the tail of the limiting scalar pdf. 

There has been a wealth of theoretical efforts addressing this last issue of the tail 
j,|l^,|ll|,||jl6|Jl|,|2^,||,|5|j2|,|2|,|^,|^ A somewhat common theme, particularly 



in the pumped case, is the prediction that the scalar pdf should develop an exponential tail. For example 
Kraichnan |^^, Shraiman and Siggia and Balkovsky and Falkovich |^ all find exponential tails. An- 
other important question is to understand the pdf of the scalar gradient. Naturally, gradient information 
may be expected to amplify contributions from small scales, and a general theory relating the scalar 
tail with the gradient tail, even for passively evolving quantities would be quite valuable. There has 
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been somewhat less theoretical effort aimed at exploring the difference in statistics between the scalar 
and the scalar gradient. Chertkov, Falkovich and Kolokolov Chertkov Kolokolov and Vergassola 
and Balkovsky and Falkovich have explored this question and have found a stretched exponential 
distribution of the scalar gradient in situations for which the scalar has an exponential tail. Holzer and 
Siggia |po| , pl| , and Chen and Kraichnan have observed similar phenomena numerically. 

In this paper we examine the scalar and scalar gradient pdf tail in an exactly solvable model first 
studied by Majda and McLaughlin and Majda who were able to construct explicit moment 
formulas for the moments of a passive scalar advected by a rapidly fluctuating linear shear flow in terms 
of iV-dimensional integrals. In that work, it was established that the degree of length scale separation 
between the initial scalar field and the fluid flow is inherent to the development of a broader than Gaussian 
pdf. 

Here, we explicitly calculate the tails of the pdf for this model. We begin by analyzing the expression 
derived by Majda for the large time 2A^th moment of the pdf for the random uniform shear model, 
which is given by an integral over R^. From these normalized moments, we will construct the tail of 
the associated pdf. We point out that in this calculation the convergence of the pdf for finite time to 
the pdf for infinite time is weak - for fixed moment number the finite time moment converges to the 
limiting moment. The convergence is almost certainly not uniform in the moments. For a more thorough 
investigation of the uniformity of this limiting process in the context of general, bounded periodic shear 
layers, see Bronski and McLaughlin 

The tail is calculated in two steps. First, using direct calculation and gamma function identities we are 
able to reduce the A^-dimensional integral to a single integral of Laplace type, from which the asymptotic 
behavior of the 2A^th moment follows easily. The asymptotic behavior of the moments is important for 
determining the tails of the probability distribution function, as we establish below. Second, we consider 
the problem of reconstructing the probability measure from the moments. Using ideas from complex 
analysis, mainly some basic facts about entire functions of finite order and type, we are able to provide 
rigorous estimates for the rate of decay of the tails of the measure. We find that the tails decay like 

exp(-c„|T|3fe) 

so depending on the precise value of the parameter a (defined in section II, below, which sets the degree of 
scale separation between the scalar and flow field) the model admits tails which are Gaussian, exponential, 
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or stretched exponential. We also show that in this model higher order derivatives of the scalar in the 
shear direction are always more intermittent, with a very simple relationship between the exponents of 
the scalar and its derivative. The distributions of derivatives in the cross-shear direction, however, display 
the same tails as the scalar itself. 

We remark that, while the stream- line topology for shear profiles is admittedly much simpler than 
that in fully developed turbulence, the fact that the exact limiting tail for the decaying scalar field may 
be explicitly and rigorously constructed suggests such models to be exceptionally attractive for testing 
the validity of different perturbation schemes. It is also extremely interesting because it demonstrates 
that, at least for unbounded flows, a positive Lyapunov exponent (as would typically occur for a general 
Batchelor flow) is not necessary for intermittency. For an interesting discussion of the role of Lyapunov 
exponents in producing intermittency see the work of Chertkov, Falkovich, Kolokolov and Lebedev. [O 



A. The random shear model 

Here, we briefly review the framework of the random shear model |2^,^,^,||. We follow Majda, and 
consider the free evolution of a passive scalar field in the presence of a rapidly fluctuating shear profile: 

dT dT 

_ + ^(i),(,)_^^Ar. (1) 

The random function, 7(i), represents multiplicative, mean zero Gaussian white noise, delta correlated 
in time: 

m^{s))=5{\t-s\) 

where the brackets, (•), denote the ensemble average over the statistics of ^{t). The original model 
considered by Majda involved the case of a uniform shear layer, v{x) = x, which leads to the moments 
considered below |^^. It a quite general fact, not special to shear profiles, that a closed evolution equation 
for the arbitrary N-point correlator is available for the special case of rapidly fluctuating Gaussian noise, 
see work of Majda |^ for a path integral representation of this fact for the special case of random shear 
layers. For the scalar evolving in (P, the N point correlator, defined as: 

^N{^.y.t)^(^T{x„y,,t^ (2) 

X = {xi,X2,X'i, ...,Xn) 



y = (2/1, y2, 2/3, ■••,yAf) 

is a function: tp^ ■ x [0, 00) satisfying 

dij^N _ . , , 1 ^ \ ^ \ ^^^^ fi^ 

= kA^^v^^ + _^;^^,(^,),(^^.)^ (3) 

where ^2N denotes the 2N dimensional Laplacian. 

We next describe the initial scalar field. Following Majda we assume that the scalar is initially a 
mean zero, Gaussian random function depending only upon the variable, y: 

T\t=o^ [ dW{k)e^^'''y\k\^Mk) a>-l (4) 

Here, (j)o{k) denotes a rapidly decaying (large k) cut-off function satisfying 4>o{k) = (p{—k), 0o(O) 7^ and 
dW denotes complex Gaussian white noise with 

(dW)^ = 
(dW{k)dW{ri))^ = 5{k + ri)dkdr] 

The spectral parameter, a appearing in (Q) is introduced to adjust the excited length scales of the 
initial scalar field, with increasing a corresponding to initial data varying on smaller scales. We remark 
that the more general case involving initial data depending upon both x and y, and data possessing both 
mean and fluctuating components, was analyzed McLaughlin and Majda 

For this case involving shear flows, the evolution of this N point correlator may be immediately con- 
verted to parabolic quantum mechanics through partial Fourier transformation in the y variable. For the 
particular initial data presented in (^), this yields the following solution formula: 

N 



^N = I e2-*-yV;^(x,k,i)TT^(fcj)|fc,|fdW^(fc,) 



where the N-body wavefunction, '0Ar(x, k, i) satisfies the following parabolic Schrodinger equation: 



dt 

V'Ar|t=0 = 1 



kAx - Vi„t(k, x)V'Ar (5) 



The interaction potential, Vmt (k, x) , is 



N 



V^nt = 4^'|k|2 + 27r2(^ kjv{xj)f 



For the special case of a uniform, linear shear profile, with v{x) — x, the quantum mechanics problem 
in (^) is exactly solvable in any spatial dimension. Taking the ensemble average over the initial Gaussian 
random measure using a standard cluster expansion, the general solution formula for {'ipN{:ic,y,t))-^^ 



is obtained |29|J3^ in terms of N dimensional integrals. The normalized, long time flatness factors, 
A*2JV = linit^oo |j2^iv 1 are calculated by evaluating the correlator along the diagonal, 

y = {v,y-,y,-- ■ -.y) 

and utilizing the explicit long time asymptotics available through Mehler's formula. This leads to the 
following set of normalized moments for the decaying scalar field, T: 

„ _ (2iv)! [ ,,. nf^ii%r 

dk- 1^1" 



/ /fi -^/cosh |fc| 

Observe that these normalized moments depend upon the parameter a. By varying this parameter 
Majda and McLaughlin established that the degree of scale separation between the initial scalar and flow 
field is important in the development of a broader than Gaussian pdf . They demonstrated through 

numerical quadrature of these integrals for low order moments that as the initial scalar field develops 
an infrared divergence (with a — > — 1, corresponding to the loss of scale separation between the initial 
scalar field, and the infinitely correlated linear shear profile) the limiting single point scalar distribution 
has Gaussian moments |29|. Conversely they showed that as the length scale of the initial scalar field 
is reduced, corresponding to increasing values of a, the limiting distribution shows growing moments 



indicative of a broad tailed distribution |32|. On the basis of these low order moment comparisons, these 
studies suggest that within these models, the limiting pdf should be dependent upon the scale separation 
between the scalar and fiow field. A fundamental issue concerns whether and how this scale dependence is 
manifest in the pdf tail. Below, we address precisely this issue, and rigorously establish that the intuition 
put forth by Majda and McLaughlin is correct through the explicit calculation of the limiting pdf tail. 



II. ASYMPTOTICS OF THE PROBABILITY DISTRIBUTION 
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A. Notation 



Recall from the previous section that the work of Majda derived exact expressions for the moments 
of a one parameter family of models indexed by the exponent a. In the remainder of the paper dfi"{T) 
will denote the probability measure for the passive scalar T in the Majda model with exponent a. The 
i*'* moment of the probability measure dfx"{T) will be denoted by fif. In this particular model the 
distribution is symmetric and thus all odd moments vanish. 



B. Large N sisymptotics of the moments 

In this model the exact expression for the 2Afth moment is given by 

„ {2Ny. /- nf^ilfcil" ,, „ 

^^2N = NnNAn / / dfeidfezdfcs ■ ■ ■ dkN 

^ Vcosh(|fc|) 



-\/cosh(fc) 

As noted by Majda cosh(|fc|) < ]^cosh(A;i) which implies the normalized flatness factors are strictly 
larger than those of a Gaussian, implying broad tails. The simplest way to analyze this integral, and 
in particular to understand the behavior for large N , is to introduce spherical coordinates. Spherical 
coordinates in N dimensions can easily be constructed iteratively in terms of spherical coordinates in A''— 1 
dimensions as follows. The coordinates in N dimensional spherical coordinates are {r, 61,62, 63. . . 6n-i}- 
If {xi~^,X2~^ ■ ■ -x^Zi} are coordinates on R^~^ then coordinates on are given by 

xf = xf-^ sin(6'Ar-i) j G 1 . . . iV - 1 
x% = r cos^jv-i 

Using this construction it is simple to calculate that the volume element in dimensional spherical 
coordinates is given by 

JV-l 

dxidx2 ...dxN = r^~^dr sm^-\ej)d6j 61 e [0, 27r] 0i>i G [0, n]. 

Since the volume element is a product measure the A'' dimensional integral factors as a product of N one 
dimensional integrals and we are left with the expression 
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(2NV 



where the Ij are given by 



7o(Ar)= /~r^(«+i)-i- 



■\/ cosh(r) 

/.27r 

h= |sin(6i)|"|cos(6')|"d6i 
Jo 

Ij = r I sin(^)p("+i)-i| cos(^)rd^ j > 1. (7) 
Jo 

The angular integrals can be done explicitly in terms of gamma functions, using the beta function 
identity 



2 r |sin(^)|2^-^|cos(0)|2»-id0 = = 
io 



which leads to the expression 



(27v)!(r(^))A^-i y r(j^) 



Observe that the product telescopes - the numerator of one term is the denominator of the next - giving 
the final expression 

(27V)! (r(^))^ fMi)-r dr 
M2„ ^^N2NNi r(Ars±i) J ^/3^ih(0 
(2Ar)! (r(H±l))A^ 
-^a^2^7V! r(iV^^±i) 

The integral over the radial variable Io{N) cannot be done explicitly, but the large A'' asymptotics are 
given by 

Io{N) « 2^(«+i)+5r(7V(a + 1)), 

so that the large N behavior of the moments is given by 

„ .,g^.+i(2iV)! r(jV(a + l))(r(^))^ 

'^^iv-^ r(7V(«±i)) • ^^^> 

Note that since 
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r(iv(a + i)) 

_____oc) as AT -.00 

the moments are strictly larger than the moments of the Gaussian. We will use this to provide rigorous 
quantitative estimates for the tails of the distribution. 



C. The Hamburger moment problem 



Having computed simple expressions for the moments of the pdf, as well as asymptotic expressions for 
large moment number, it is natural to ask the question of whether one can do the inverse problem, and 
deduce the pdf itself. The problem of determining a measure from its moments is a classical one, known 
as the Hamburger moment problem |^,^,^. This problem has a rich theory, and we mention only a 
very few of the most basic results here. For an overview of the subject, see the book by Shohat and 



Tamarkin or the recent electronic preprint by Simon |44 



The two most important questions are, of course, existence and uniqueness. There is a necessary and 
sufficient condition for a set of numbers {fJ-i} to be the moments of some probability measure, namely that 
the expectation of any positive polynomial be positive. This translates into the following linear algebraic 
conditions on the diagonal determinants of the Hankel matrix, the matrix with entry Hi+j: 



|a*o| > 0, 



m fJ-2 



> 0, 



> 0. 



Mo /^l M2 
Ml M2 M3 
M2 A*3 M4 

These conditions appear to be quite difficult to check in practice. However since the moments considered 
here are, by construction, the moments of a pdf this condition must hold. 

A more subtle question is the issue of uniqueness of the measure, usually called determinacy in the 
literature of the moment problem. One classical sufficient condition for the determinacy of the moment 
problem is the following condition, due to Carleman : If the moments /Lt„ are such that the following 
sum diverges 



then the Hamburger moment problem is determinate. Given the asymptotic expression for the moments 
given in Equation ( |ic| ) it is easy to check that 
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and thus there is a unique measure with these moments for — 1 < a < 1. We wiU see later that this 
corresponds to probabihty distribution functions with tails that range from Gaussian through exponential. 

In the case a > 1 which, as we will see later, corresponds to stretched exponential tails, the problem 
probably does not have a unique solution. Indeed there are classical examples of collections of moments 
with the same asymptotic behavior as the stretched exponential distribution for which the moment 
problem has a whole family of solutions. 

Given this we come to the question of actually calculating the measure given the moments. There 
is a rather involved theory for this in the determinate case involving, among other things, orthogonal 
polynomials and continued fractions p7|p^ , but in general this problem is extremely difficult. However 
we show in the next section that it is relatively straightforward to reconstruct the tails of the measure 
from the moments. 



D. Asymptotics of the tails of the distribution 

Recall that the 2A^th moment of some probability measure dfj,°'{T), 

f^^N = I T'^'dfi-iT). (11) 

We are interested in calculating the asymptotic rate of decay of the tails of the probability measure dfi" (T) . 
The information about the behavior of the tails of the distribution is contained in the asymptotic behavior 
of the large moments. We study the tails of the measure dfi"{T) by introducing the function 

where C is some as yet unspecified constant. The factor of r( ) jg chosen so that the series for /" 

has a finite but non-zero radius of convergence. This will give us the sharpest control over the tails of 
dfi°'{T). It is convenient to demand that the radius of convergence of the series be one. Using the root 
test it is straightforward to check that the radius of convergence of the sum is given by 

* _ po-(a+2) (« + 3)^ / ^ 

(a + i)^Vrm' 

so we choose 



10 



(a + 3)^V ^ 

Since the coefficients /i2jv s-'^s tfie moments of a probability measure d(i°'{T) we tiave tfie alternative 
expression 

When z is inside the radius of convergence of the sum (i.e. \z\ < 1) we can switch the integration and 
the summation and get the following expression for /" 

rw-/gj3;p^*°m (13) 

= J F°'{zT)4i°'{T). (14) 

We note a few simple facts. First notice that the function f"{z) is a kind of generalized Laplace transform 
of the measure dfj,°'{T). The quantity inside the integral, F°'{zT) = ^ (g2jr(^J<^^+°> ) *^onverges absolutely 
for all z and thus F"{zT) is an entire function of the complex variable z. Further we know that the 
integral must converge for |z| < 1 and diverge for some \z\ > 1, since the original series converged in a 
circle of unit radius. We note that the entire function F"(z) satisfies 

<F-{\^\). (17) 



so the function F"(z) grows fastest along the real axis. Thus we know that the integral in Equation (14) 
converges for — 1 < z < 1 and diverges for z > 1, z < —1. Thus the problem of understanding the rate of 
decay of the tails of the probability measure 4^"(T) has been reduced to that of determining the rate of 
growth of the function F{zt). There is a well-developed theory for studying the rate of growth of entire 
functions, the theory of entire functions of finite order. We recall only the basic facts here - the interested 
reader is referred to the texts of Ahlfors ^ and Rubel with Colliander . 

The radial maximal function Mp{r) of an entire function F{z) is defined to be the maximum of the 
absolute value of F over a ball of radius r centered on the origin: 

Mpir) = max |F(z)| 
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The order p of a function F is defined to be 

,. log+log+MF(r) 
p = hmsup , 

where log^(a;) = max(0, log(x)), if this limit exists. It is easy to see from this definition that F is of 
order p means that F grows asymptotically like exp(A(z)|z|'') along the direction of maximum growth in 
the complex plane, where A(z) grows more slowly than any power of z. A related notion is the type of a 
function of finite order. If F is of order p then the type r is defined to be 



r = limsup 



rP 

when this limit exists. Again speaking very roughly the type r gives the next order asymptotics: if F is 
of order p and type r then F grows like B{z) exp(T|z|''), where B{z) is subdominant to the exponential 
term. Note that by Equation (|l^) the function F" grows fastest along the real axis, and thus the maximal 
rate of growth in the complex plane is exactly the rate of growth along the real axis. 

There exist alternate characterizations of the order and type of a function in terms of the Taylor 
coefficients An which are very useful for our purposes. These are given as follows: 

log+ \og+ Mpjr) nlogjn) 



hmsup — = — limsupri|ii„|^' . 



T = nmsup '■ = — nmsupri|A„|''' '". (19) 

r — *oo pG n — 'OO 

For the proofs we refer to the text of Rubel with CoUiander . Using the expressions given in equations 

18|) and ( [l^ ) we find that the order p and type r of F"(2) are given by 

„ 2nlog(2n) 4 

p = hmsup — 



log(C2»r(i2^)) 3 + a 

1 , , (3 + a)n , , -p 1 

r" = -limsupn|r(^^)|- = - 

pC n— too ^ 



Thus we know that F"{zT) grows like A{zT)exp{C ''|z|3+° |r|3+=) along the real axis, where A{zT) 

4 

grows more slowly than exp(Z?|r|3+° ) for any D. Further we know that the integral 

J F°'{zT)dp°'{T) 

converges for \z\ < 1 and diverges for z > 1 or z < — 1, so to leading order the rate of decay of the measure 
dpf^iT) is given by exp(— |C|^'*/^'^+"^|r|'*/*^'^+"''). It is easy to check that as a — + — 1 this estimate becomes 
exp(— recovering the normalized Gaussian. 
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This result is probably best restated in terms of the cumulative distribution function, rather than the 
probability measure. If P{T, T') = dn{T), with T' > T, then it is easy to show that the above implies 
that 

lim ex.p{c\T\^)P{T,T') =0 c<\C\^ 

T— >oo 

= 00 c> \C\'^ 



III. INTERPRETATION AND CONCLUDING REMARKS 

Physically the Majda model can be thought of as a model for the behavior of a passive scalar at small 
scales, when the scale of the flow field is much larger than the scale of the variations of the scalar. Recall 
that the random scalar is given by 

T{y) = j \k\^Mk)e^"'^'>dW{k) (20) 
< T{y)T{y') > =J \knMk)\^e^^"'^y-y'Uk. (21) 

In the limit as a approaches —1 there is an infrared divergence, so that the energy of the scalar is 
concentrated at larger and larger scales. In this case 2, so the normalized distribution function 

becomes Gaussian, as was originally observed by Majda. 

One important fact about this model which we would like to emphasize is that it predicts that higher 
derivatives of the advected scalar should be increasingly intermittent, a fact which was not strongly 
emphasized in previous work. Observe that due to the special nature of shear flows the scalar derivative 
dT/dy satisfles the same equation as the scalar T with no additional terms!. We further note that the 
initial condition for the derivative of the scalar is given by 

8T f 

-g:^= J 27ri\k\'^k(l)o{k)e'^'''''^dW{k) (22) 
< > =47r'J \kr+'\Mk)\'e'''"'^y-y'Uk, (23) 

so the derivative of the scalar has a representation of the same form as the representation of the scalar 
itself, but with the exponent a increased by two, and a slightly modified (po{k). Recall that the exponent 
a determines the amount of energy at the largest scales and thus the degree of intermittency, with the tails 
decaying as exp(— r^/(^+")). Our calculation shows that increasing a increases the width of the tails of 
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the probability distribution function, implying that derivatives are more intermittent! These predictions 
for the behavior of the tails of the scalar as compared with the scalar gradient are in extremely good 
agreement with experimental and numerical results. For instance our calculation shows that if the scalar 
has exponent a = — 1, so that the probability distribution function of the scalar has Gaussian tails, then 
the derivative of the scalar has exponent a = 1, implying that the distribution of the derivative has 
exponential tails. This agrees quite well with the experiments of Van Atta and Thorddsen jsot , as just 
one example, who observe that in turbulent thermally stratified flow that the pdf for the density has 
Gaussian tails, while the pdf of the density gradient has exponential tails. Similarly if the scalar has 
exponent a — 1, so that the distribution of the scalar itself is exponential, then derivative of the scalar 
should have exponent |. This agrees with the calculations of Chertkov, Falkovich and Kolokolov p^ , 
and Balkovsky and Falkovich |^ also predict exponential tails for the scalar and stretched exponential 
tails with exponent | for the scalar gradient in the Batchclor regime. This also shows reasonably good 
agreement with the numerical experiments of Holzcr and Siggia po| , pl| . In their experiments Holzer and 
Siggia find that a scalar with exponential tails has a gradient with stretched exponential tails. For large 
Peclet number the exponent of these stretched exponential tails is in the range of .661 — .563. 

Of course one can eliminate a entirely, and one finds the following relationship between the distribution 
of the scalar and the scalar gradient within this model. If T is distributed according to a stretched 
exponential pdf with exponent p, and the gradient Ty according to a stretched exponential pdf with 
exponent p', then p, p' are related by 

1 11 

2 + p " 7' 

It would be extremely interesting to check if this relationship, or some generalization of it, holds in greater 
generality than shear flows. The above numerical and experimental evidence suggest that this might not 
be an unreasonable hope. 

The distribution of the x, or cross-shear, derivatives can also be calculated using the same explicit 
representations derived by Majda. Calculations by the authors for deterministic initial data have shown 
that derivatives in the cross-shear direction have a distribution with the same asymptotic behavior as 
the scalar itself. This should be compared to and contrasted with the papers of Son , and Balkovsky 
and Fouxon Q, which predict distributions with very broad tails (all of the higher moments diverge as 
t — > oo) and which predicts the same distribution for derivatives of the scalar as for the scalar itself. 
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We would also like to comment on the relationship between intermittency and the Lyapunov exponents 
of the underlying flow field. A number of papers have addressed the problem of intermittency in the large 
Peclet number limit by attempting to relate broader than Gaussian tails to the Lyapnuov exponents of 
the flow field jl^ . It is worth noting that a shear flow does not possess a positive Lyapunov exponent, 
but as we have shown here a shear flow can generate exponential and stretched exponential tails in the 
passive scalar. This shows that chaotic behavior in the underlying flow, while probably an important 
effect in realistic flows, is not necessary for the generation of broad tails and intermittency. 

Finally we would like to comment on the rate of approach to the limiting measure in time. The results 
presented here analyze the infinite time limit of the measure. As mentioned earlier the convergence to 
this limiting measure is expected to be highly non-uniform. A preliminary calculation by the authors for 
a special choice of the cut-off function (/)o(fc) suggests that for large but finite times the pdf looks like the 
pdf for the infinite time problem in some core region, with Gaussian tails outside this core region. As 
time increases the size of this core region demonstrating non-Gaussian statistics grows, and the Gaussian 
tails get pushed out to infinity. We believe this same picture to hold for any choice of the cut-off function 
<j>o{k), but more work is needed to establish this fact. 
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